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ABSTRACT 

Type 2 quasars are luminous Active Galactic Nuclei whose central engines are seen through large 
amounts of gas and dust. We present Spitzer spectra of twelve type 2 quasars selected on the basis of 
their optical emission line properties. Within this sample, we find a surprising diversity of spectra, from 
those that are featureless to those showing strong PAH emission, deep silicate absorption at 10/im, 
hydrocarbon absorption, high-ionization emission lines and H2 rotational emission lines. About half of 
the objects in the sample are likely Compton-thick, including the two with the deepest Si absorption. 
The median star-formation luminosity of the objects in our sample measured from the strength of 
the PAH features is 5 x lO^^L©, much higher than for field galaxies or for any other AGN sample, 
but similar to other samples of type 2 quasars. This suggests an evolutionary link between obscured 
quasars and peak star formation activity in the host galaxy. Despite the high level of star formation, 
the bolometric output is dominated by the quasar in all cases. For a given strength of 10/im Si 
absorption, ULIRGs are significantly colder than are type 2 quasars (their [14.5/im]/Fiy[27.5/im] 
ratio is 0.5 dex lower), perhaps reflecting different obscuration geometries in these sources. We find 
that the appearance of the 10/im feature (i.e., whether it shows in emission or in absorption) is well- 
correlated with the optical classification in type 1 and type 2 quasars, contrary to some models of 
clumpy obscuration. Furthermore, this correlation is significantly stronger in quasars (Lboi ^ 10*^ 
erg/s) than it is in Seyfert galaxies (iboi ^ 10^^ erg/s). 

Subject headings: galaxies: active - galaxies: nuclei - infrared: galaxies - quasars: general 



1. INTRODUCTION 

Type 2 quasars are powerful (Lboi ^ 10"*^ erg/s) ac- 
tive galactic nuclei (AGNs) whose central regions are 
shielded from the observer by large amounts of gas and 
dust. Until recently, the existence of dust-enshrouded 
quasars was in question, but now these objects 
are r o utinely discovered at X-ray fe .g.. ISzokoly et aTl 
20041 iNetzer et al.' l2006l), inf rared (iLac v et all l200l 



Stern et al. 2005; Martfncz-San sigre et al< 2006) and op- 



tical (Zakamska ct al. 2003j; iReves et al.l 12008)) wave- 
lengths. A significant controversy remains regarding 
their space density and contribution to the AGN pop- 
ulation luminosity budget. Some analyses suggest that 
there are at least as many obscured as unobscured 
quasars even at the highest luminosity end (Reyes et alj 
[2OO8; Martmez-Sansigre et al. 2006; Polletta^eXal,._20Qj), 
whereas others argue for a strong decline of ob- 
scured quasar fraction with luminosity (,Ueda et al.|[2003t 
iMaiolino et"alll2007t P&eister et al.ll2008[ K 

If an AGN is entirely embedded in an optically thick 
dusty cloud, virtually all of its original radiation (emit- 
ted primarily at optical to soft X-ray wavelengths) is 
absorbed by dust and is then re-emitted thermally in 
the IR. In this case, there is no optical signature of the 
AGN, but its presence may be revealed when hard X- 
rays and a large IR luminosity are detected. In this way, 
some ultra-luminous infrared galaxies (ULIRGs ) have 
been shown to contain embedded AGNs (.Risaliti et al.l 
I2OOOI: iFranceschini et all 120031 iTeng et all l2005l ) . The 
dominant energy source in ULIRGs remains a matter 



of controversy. Some authors argue that they are pow- 
ered by star form ation, with the AGN contributing only 
a small fraction (jFranceschini et a l. 2003), while others 
suggest that bu ried AGNs are energ etically dominant in 
many ULIRGs (jlmanishi et al.ll2007l) . 

An alternative possibility is that the nucleus is only 
partially covered by the obscuring material. In this case, 
the AGN is seen as obscured or unobscured depending 
on the line of sight, making its UV, optical and X-ray 
emission highly anisotropic (jAntonuccil Il993l ) . In par- 
ticular, UV and X-ray emission can make its way out 
through the holes in the obscuring cloud and photo- 
ionize extended regions, which in turn produce power- 
ful optical emission lines. Using the Sloan Digital Sky 
Survey (SDSS: lYork et al l l200a [Sto ughton et al.l l200l 
lAdelman-McCarthv et al.l l2008l ) we have selected hun- 
dreds of type 2 quasars at redshifts z < 0.8 on th e basis 
of their optical emission line properties ( Zakamska et al.l 
[200I IReves erall[2008h . 

IR observations are key to determining the physical 
conditions of the obscuring material and its geometric 
properties such as its covering factor and extent. In this 
pape r we present Spitzer Space T elescope ( Werne r et al.l 
I2OOI Infrared Spectrograph fIRS: lHouck et al.il200l ob- 
servations of twelve type 2 quasars from the SDSS sam- 
ple. In Section [2] we describe our sample selection and 
data reduction. In Section [3] we present a detailed analy- 
sis of the spectral features. We compare our sample with 
other samples of type 2 quasars in Section |4l we discuss 
our results in the context of radiative transfer models in 
Section [5] and we conclude in Section [6l 
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Throughout the paper, we adopt a cosmology with 
h — 0.7, Qm — 0.3, r^A — 0.7. Objects are identified 
as SDSS Jhhmmss.ss+ddmmss.s in Figure [T] and Table [T] 
and are shortened to SDSS Jhhmm+ddmm in the text. 
We refer to type 1 and type 2 AGNs based on their op- 
tical classification (in type 1 objects, broad lines are de- 
tected and in type 2, they are not). When discussing the 
luminosity dependence of AGN phenomenon, we distin- 
guish between Seyfert galaxies (Lboi 10^^ erg/s) and 
quasars (Lboi ^ 10^^ erg/s) based on their bolometric 
luminosities. 

2. SAMPLE SELECTION, OBSERVATIONS AND DATA 
REDUCTION 

The type 2 quasars discussed in this paper were se- 
lected from the spectroscopic database of the SDSS 
on the basis of t heir op tical emission l ine proper- 
ties ( Zakamska et a l. 2003; Hao et al.ll2005aHReves et all 
|20G8). Briefly, type 2 AGNs are required to show narrow 
emission lines with ratios characteristic of a non-stellar 
ionizing continuum. We then use the luminosity of the 
[OIII]5007A emission line as a pr oxy for the intrinsic lu- 
minosity of the obscured nucleus (jZakamska et al.ll2003f) 
and classify type 2 AGNs as quasars if L[OIII]> 10^-^ 
and as Seyfert galaxies otherwise. iReves et al.l ()2008f ) 
obtained a sample of 887 type 2 quasars at redshifts 
< z < 0.8 by applying these criteria to the most recent 
SDSS data. [OIII] luminosities were not corrected for 
reddening, and they are thus likely to be underestimated 
(IRe ves et al . 2008). 

Our extensive follow-up work has demonstrated that 
the objects that we classified as type 2 quasars based 
on their optical properties indeed fit operational defini- 
tions of type 2 quasars used at other wavelengths. In 
particular, they typically have high X-ray luminosities, 
> lO**^ erg/s, and large amoun ts of neutral gas 
along the line of sight, Nh > 10^^ cm~^ (IZakamska et al.l 
[2001 iPtak eral]l2006l : iVignah et all 12004 120061 ^. Fur- 
thermore, using optical spectropolarimetry, we demon- 
strated that at least some of the objects in our sam- 
ple harbor broad-line (type 1) AGNs in their centers 
(jZakamska et al.ll2005l ). 

In Cycle 1 of the Spitzer Space Telescope we con- 
ducted seven-band IRAC+ MIPS photometry of 25 type 
2 quas ars from IZ.akamska et al.l (|2003f l and iHao et all 
(|2005af) . These objects were selected primarily based on 
their [OIII]5007A line luminosity (specifically, L[OIII]> 
IO^Lq). Spitzer photometric observations (Zakam- 
ska et al., in prep.), as well as archival IRAS data 
(IZa kamska et al.ll2004[ ). show that type 2 quasars from 
our sample typically have IR luminosities in excess of 
10^^ erg/s, placing them into the classical quasar lumi- 
nosity regime and validating our [OIII]5007A-based lumi- 
nosity criterion. We selected ten sources from this sam- 
ple with the highest IRAC-8^m and MIPS-24^m fluxes 
(i^j/[8/xm] > 1.5 mJy and F,y[24^m] > 6 mJy) for spec- 
troscopic follow-up with the IRS (Table [T]). 

In addition, we matc hed the most rece nt catalog of 
SDSS type 2 quasars (jReves et al.l l2008f ) against the 
Spitzer archive and found the following two sources with 
archival IRS observations: 

1. SDSS J091345.49-h405628.2 (=IRAS 

B09104-f4109) - program 1018 (In-orbit Checkout 
/ Science Verification), PI Armus, z=0.442; 



2. SDSS J005009.81-003900.6 (=IRAS B00476-0054) 
- program 105 (Guaranteed Time Observations), 
PI Houck, z=0.729. 

The first of t hese sources is a wel l-known luminous 
type 2 quasar (jKleinmann et allll988[ ). whereas the sec- 
ond source was originally classified as an ultra-luminous 
infrared galaxy (ULIRG). Both objects show high- 
ionization optical line ratios, which is why t hey were se- 
lected by the optical procedure described in IReves et al.l 
(2008), so we include them in our analysis. 

The ten sources in our targeted program were observed 
with the first order of the short wavelength low reso- 
lution module (SLl; 484 sec exposure for each object) 
and with the second order of the long wavelength low 
resolution module (LL2; 488 sec). These observations 
cover the wavelength range 7.4/^m < Aobs < 21.3/im 
with resolution i? ~ 57 — 127. We started our analy- 
sis from the Basic Calibrated Data (BCD) provided by 
the Spiteer Science Center (pipeline ver. S13.2.0), which 
applies dark corrections, flat fields, and droop and non- 
linearity corrections, and provides flux calibration. We 
removed the sky background and most of the hot pixels 
by a pairwise subtraction of frames at different nod po- 
sitions. We then extracted positive and negative spectra 
from the subtracted images with the SMART package 
(iHiedon et al.ll200l . using the 'interactive/column/No 
SkySub' method. Additionally, we applied hot pixel cor- 
rections to the subtracted two dimensional spectra us- 
ing IRSCLEAN^. The resulting spectra (4 for SLl and 
8 for LL2 for each object) were averaged into the final 
spectra shown in Figure [TJ When necessary, the spec- 
tral errors were estimated from pixel-to-pixel variations 
by subtracting a smooth model from the spectrum. We 
found that our reduced spectra are in agreement with 
IRAC-8.0 and MIPS-24 photometry and conclude that 
the absolute flux calibration is accurate to better than 
15%. 

We reduced the two archival sources in the same way 
(pipeline ver. S14.0.0), and they are also shown in Fig- 
ure [TJ The archival observations had total exposures 2-4 
times shorter than the targeted observations, since the 
archival targets are IRAS-detected and therefore are rel- 
atively bright in the IR. We have included SLl, LL2 and 
LLl observations (total wavelength coverage 7.4/im < 
Xobs < 38.0^m). In the source SDSS J0913-I-4056 
(=IRAS B09104+4109), one of the two available LL2 
datasets produced fluxes which were systematically lower 
(by about 5%) than those from the other LL2 dataset 
and those extrapolated from the LLl dataset. We multi- 
plied the fluxes of the discrepant LL2 dataset by a linear 
function of wavelength to correct for this difference. 

Some of the analysis in this paper makes use of IR 
photometric data. For the ten targeted sources, the 
IR photometry is from our Cycle 1 Spitzer program 
(Zakamska et al, ,in prep.). For SDSS J0913-f4056 
(=IRAS B09104+4109), we use IRAS catalog fluxes^ 
(^^^[12^m] = 130 mJy, F^[25fj,m] = 333 mJy, 
i^y[60/im] — 525 mJy) and Spitzer fluxes derived from 
IRAC observations in program 105. Using aperture pho- 
tometry, we calculate F^[3.6fini] = 5.0 ± 0.7 mJy. For 

^ http: / /ssc. spitzer. caltech.edu/archanaly/contributed/irsclean 
taken from the NA SA/IPAC Infrared Science Archive, 
|http://irsa.ipac.caltech.edu[ 
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SDSS J00 50-0039 (=IRAS B004 76-0054), we use IRAS 
data from lStanford et all (|2000D (i^^[60^m] = 260 mJy). 
In all cases where the data overlap, the IR photometry is 
in agreement with the fluxes derived from the IRS spec- 
tra within the uncertainties of the photometric data. 

3. ANALYSIS OF SPECTRAL FEATURES 

3.1. Identification and detection of spectral features 

From Figure[T]it is apparent that the IR spectra of type 
2 quasars display a wide range of properties. Some are 
nearly featureless, whereas others are busy with absorp- 
tion and emission features, some of which are identified 
in the top panels of Figure [TJ The difficulties of placing 
the continuum to identify detected features, as well as of 
disentangling blends of em ission and absorption arc well 
known for IR spectra (e.g.. lSpoon et al. 2002. 2QQA), and 
some of our conclusions may suffer from the uncertainties 
in this rather subjective step. Nevertheless we attempt 
to minimize these uncertainties or at least obtain their 
measure by using several different approaches to spectral 
analysis. 

In the first step, we use PAHFIT (jSmith et al.ll2007f) . a 
public IDL-based routine that decomposes low-resolution 
IRS spectra into a combination of starlight (represented 
by a black body spectrum), featureless continuum (rep- 
resented by a sum of thermal components at several 
fixed temperatures), emission features (PAH dust emis- 
sion features, fine-structure lines of ions, rotational lines 
of II2), and dust extinction (modeled as a screen of cold 
silicate dust affecting all other components). The code 
was developed primarily for the analysis of polycyclic 
aromatic hydrocarbon (PAH) features in high signal-to- 
noise ratio spectra of star forming galaxies, so the rela- 
tive strengths of PAHs are allowed to vary. The advan- 
tage of using this code for our objects is that it produces 
excellent global fits in every case. The disadvantage is 
that the large number of fit parameters (30-40) makes it 
hard to assess which of the weak features included in the 
best PAHFIT model are actually detected. Even for the 
stronger features, s ome fit p arameters are degenerate, as 
was illustrated by Spoon et al. (2002) in the case of the 
interplay between PAH emission and silicate (Si) absorp- 
tion. We use PAHFIT as a guide to indicate significant 
PAH feature detections. 

In the second step, we use a much simpler global fit 
to the data consisting of a featureless continuum (repre- 
sented by a second order polynomial) and a PAH tem- 
plate with fixed ratios of feature strengths. The sum 
of these components is then multiplied by exp[— t(A)], 
where t(A) is the silicate opacity curve, whose normal- 
ization is a free parameter which is allowed to be positive 
or negative. The PAH template that we use is derived us- 
ing PAHFIT from a composite spectrum of galaxies from 
Smith et al. (2007; sohd line in FigurelU top), and the Si 
opacity curve is ba sed on ISO o bservations of the Galac- 
tic center (Kempe r et al.l [20041 : dotted line in Figure [U 
top). In some cases we masked the strongest emission 
hues in the fitting (e.g., [NeVI]7.63^m, H2 S(3)9.67^m, 
[SIV]10.5lAim and [NeII]12.8^m in SDSS J1641-h3858). 
This model has only five free parameters, but still pro- 
duces good global fits for our spectra. 

Using these two models, we calculate the Si feature 
strength, 6*9.7 = ln(/;y[9.7/^m]//jy[cont]), where fi,[9.7fiin\ 
is the model flux density at rest-frame 9.7/im and f^, [cont] 



is the flux density extrapolated from the model contin- 
uum outside the feature. In Table [T] and Figure [U we 
quote Si strengths from the second method. The la un- 
certainties shown in Figure [2] are taken to be half the 
difference in Si strengths measured by the two methods. 

We tried a range of model Milky Way and 
Small Magellanic Clo u d (SM C) opacity curves from 
IWeingartner fc Draini (1200 iD. as well as the Sgr A* 
opacitv curve from lKemper et al.l (I2004D an d the Galac- 
tic dust opacity curve from lChiar fc TielensI (|2006l ). The 
variation in 5*9.7 measured using different opacity curves 
amounts to no more than 10%, although the conversion 
between IR optical depth and optical extinction or hydro- 
gen column density varies by large factors (r[9.7//m] = 1 
corresponds to Ay = 12.5 mag for the Milky Way opacity 
curve and to Ay — 4.7 mag for the SMC curve). We flnd 
that the shape of the lO/^r n Si feature in type 2 quasars 
is well-represented by the iKemper et al.l (|2004[ ) opacity 
curve, which we use hereafter. 

In the third step we examine detections of emission 
lines, such as those from fine-structure transitions in 
ions (Section 13. 5[) and those from rotational transitions 
of molecular hydrogen (Section 13. 6p . As shown in the 
to p panel of Fi g ure [Tl we consider all transitions listed 
in ISmith et all (|2007l) . as weU as [NeVI]7.63^m and 
[NeV]14.33/xm. All these lines are expected to be un- 
resolved, so their widths are determined by the instru- 
mental resolution (FWHM^ =0.1/xm for SLl, 0.14/im for 
LL2 and 0.34/zm for LLl in the observed frame). For each 
line, we calculate the error per pixel by fitting a low or- 
der polynomial to a cut-out from the spectrum excluding 
a region 3xFWHMfl/2.35 wide on both sides of the line 
centroid. We then include the region containing the line 
into the fit and examine the reduced with or without a 
fixed-width Gaussian representing the line to determine 
the confidence of the detection. Similarly, we examine de- 
tections of PAH emission features by fitting a second or- 
der polynomial -I- fi xed-width Prude ( damped harmonic 
oscillator) profiles ([Smith et all 120071) . The line fluxes 
that we obtain from these simple fitting procedures are 
then compared to PAHFIT line fluxes, allowing us to 
estimate measurement confidence and uncertainty. For 
unresolved lines, the uncertainties are in good agreement 
with those calculated from the measured error per pixel. 

3.2. Silicates 

Silicates are a major constituent of interstellar and cir- 
cumstellar dust. Their presence is indicated in many 
classes of Galactic and extragalactic objects by two 
strong features, one due to a Si-0 stretching mode 
and peaking at around 10 /im and another due to a 
0-Si-O bending mode an d peaking at around 18/im 
([Knacke fc ThomsonI [19731) . 

In a recent study of about 200 ULIRGs, Seyfert 2 
galaxies, Seyfert 1 galaxies and type 1 quasars.l Hao et all 
(2007) find that there exists a strong relationship be- 
tween the optical classification of these objects and their 
IR properties. We show their sample together with 
our type 2 quasars and a compilation of ULIRGs and 
AGNs from different datasets in a Si strength vs. mid- 
IR color diagram (Figure [21 model lines shown in this 
Figure are described in detail in Section [5]) . For exam- 
ple, the Si feature always appears strongly in absorption 
in ULIRGs, implying a deeply embedded source of en- 
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ergy with the temperature stee ply declining toward the 
observer ([Levenson et alJ l2007l ) . The total amount of 
obscuring material is such that is optically thick even 
in mid-IR. On the contrary, it was recently confirmed 
that type 1 quasars show the Si feature in emission 
(jSiebenmorgen et al.l[2005l : iHao et al.ll2005bf l. The emis- 
sion feature is thought to be produced by the thin hot 
layer of dust on the inner side of the obscuring material, 
where it is illuminated and heated by the central engine. 
The hot dust emission region is seen without intervening 
absorption in type 1 quasars. In type 2 quasars from our 
sample the silicate feature centered at 9.7/Ltm is detected 
in eight spectra in absorption, and in the remaining four 
it is not detected, with upper limits of Sg,r < 0.1. As 
shown in Figure [2l the same is true for other samples of 
type 2 quasars. 

A possible exception is SDSS J0920+4531 (Figure O. 
This object has a hump in its mid-IR spectrum cen- 
tered at about ll/itm, which may be due to Si emission. 
With this assumption, we estimate the color tempera- 
ture of this emission - that is, we fit the continuum- 
subtracted spectrum F^(X) with thermal emission mod- 
els B^{X, T)r(A), where Bi,{\, T) is the black-body spec- 
trum with temperat ure T. Our best fit (T = llOK with 
iKemper et al.|[200^ opacitv curve) is shown in Figure [3l 
This temperature is significantly lower than those typ- 
ically der ived for type 1 q uasars showing Si emission 
(~ 200K. iHao et al.l l2005bD : such a low temperature is 
required to account for the shift of the feature centroid 
to long wavelengths. The peak strength of the emission 
is ^ii/im = 0.13, but there is no flux excess at 9.7/im 
due to Si emission, so we use 59.7 = 0.0 for this object. 

We now investigate the relations hip betwe e n the 
strength of Si feature and IR colors. iHao et all ()2007l ) 
find a strong correlation between the depth of the Si fea- 
ture and the log(i^i/[14.5/im]/F,y[27.5/im]) color in their 
compilation of ULIRGs, Seyfert galaxies and type 1 
quasars, as shown in Figure [D We therefore look for 
relationships between the values of 6*9.7 and IR colors 
within our sample of type 2 quasars. To this end, for 
each object we compile all available IR photometric data 
from our iSpiizer observations (Zakamska et al., in prep.), 
from IRAS where applicable (Section [2]), and from the 
IRS spectra themselves (from which we extract fluxes 
at rest-frame 14.5/im). We then model the observed 
spectral energy distribution as a third-order polynomial 
in the log(A) — log(Fj^) space, interpolate to calculate 
model fluxes at rest-frame 3.0, 4.0, 5.5, 8.0, 14.5, 27.5 
and 50.0/xm, and compute all 21 rest-frame flux ratios. 
We then calculate the Spearman rank correlation coefh- 
cient {ts) between values of 5*9.7 and these IR colors for 
our sample of twelve objects. We find that the values 
of rg range from 0.2 to 0.75, depending on the color in 
question, with the most common value around 0.5. The 
sense of these correlations is that the Si 9.7/i m absorp- 
tion f eature is stronger in colder objects, as iHao et aTl 
([2007^ found. Most correlations are rather weak (the 
values rg — 0.5 imply that there is a < 10% chance 
to find such relationships in uncorrelated datasets), but 
two are significant at a > 98% level. These are 6*9.7 vs. 
F^[A.{)^im\/F^[%.Q^im\ and 5*9.7 vs. F^[^.b^iTa\/ F^[b0^im\; 
they are shown in the bottom panels in Figure [2] In Sec- 
tion [5] we discuss the implications of these findings and 
the geometrical and physical factors that determine the 



positions of objects in this diagram. 

It would be interesting to determine the relationship 
between gas column density and dust column density in 
the obscuring material in AGNs. The column density of 
neutral gas can be probed directly by X-ray observations. 
The strength of the Si feature would be a direct measure 
of the dust optical depth only if the dust were so cold 
that its emission at 10/im were negligible. In practice 
dust exists at a range of temperatures and both emits 
and absorbs in the mid-IR, so that 69.7 reflects some 
combi nation of o ptical depth and temperature distribu- 
tion. IShi et all (j2006f ) find a correlation (with a large 
scatter) between the strength of the Si feature and the 
neutral gas column density measured using X-ray obser- 
vations. For 11 out of 12 objects in our sample, pointed 
Chandra and XMM observations are available (Table [1]). 
In most cases, only a few photons are detected, and we 
cannot derive a reliable measurement of the column den- 
sity from X-ray spectra. Instead, in Figure [4] we show 
the strength of the Si feature vs. the X-ray flux of our 
objects, normalized to their [OIII]5007A flux presumed 
to be an isotropic measure of the bolometric luminosity. 
We calculate X-ray fluxes expected in a type 1 quasar 
based on the [OIli]5007A flux. The mean (thick solid 
line) and the la range (dotted lines) of the Lx/i[OIII] 
ratio are taken from Hcckman et al. (2005). We then use 
the observed X-ray/ [OIII]5007A ratio as a crude mea- 
sure of the neutral gas column density (the grey lines 
show the expected ratio for two values of the column den- 
sity) . We see that half of our sample li kely has very high 
column densities, Nr > lO'^^ cm^^ (jPtak et al.l l2006t 
IVignali et al.l 120061 ). including the two objects with the 
deepest Si absorption. Although our method of estimat- 
ing Nh is rather crude, we also see that type 2 quasars 
occupy a similar range of the 5*9.7-iVg parameter space 
as do Seyfert 2 galaxies fjShi et al.ll2006f ). 

3.3. Other absorption features 

Mid-IR spectra of some deeply absorbed sources, 
both Galactic (e.g., protostars) and extragalactic (e.g., 
ULIRGs), display abso rption due to organic molecules 
and inorganic ices fe.g.. lSpoon et al1 [2002 and references 
therein). As we show in Figure [3 in SDSS J0815-f4304 
we detect a broad absorption feature between 6 and 
8/im and within it a pair of narrow features - an 
absorption complex quite simil ar to the one see n in 
the ULIRGs IRAS F00183-7111 (iSpoon et al.l [200l and 
IRAS 08572-1-3915 (iDartois et al.ll2007^ . 

We identify the two narrow absorption features cen- 
tered at 6.85/im and 7.27/im with bending (deformation) 
modes of CH bonds in the CH2 and CH3 subgroups of 
the aliphatic (chain- like) component of organic molecules 
within the dust. In Figure [5^, we compare the observed 
absorption with the opacity curve of hydro genated amor- 
phou s carbon measured in the labor atory (IDartois et all 
1200*71). Using normaliza tions from ISpoon et al.l ()2004[ ) 
and lbartois et al] (|2007f ) and peak optical depths of the 
two features (r[6.85/im] = 0.19 and r[6.85^m] = 0.12), 
we find column densities of A^cHa — A^CHa — 6 x 10^^ 
cm~^. This number is quite high; indeed, taking the 
abundance of carbon to be 3 x 10"'* relative to hydro- 
gen and assuming that 10% of it is locked into aliphatic 
chains (|Dartois et al.l [20071 ) . we get a neutral hydrogen 
column density of Nh ~ 4 x 10^^ cm~^. This value 
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may be compared with Nh ^ 10^^ cm^^ estimated from 
the lack of an X-ray d etection in this object (Figure H] 
and lVignah et al]|2006[ ). The difference between the two 
measurements is indicative of large uncertainties in all 
the steps of our calculation (such as the conversion be- 
tween optical depth of the features and -/Vch2; -^CHs, 
as well as carbon abundance and the fraction of carbon 
locked in the chains) . It is interesting that the peak opti- 
cal depth at 6.85/im that we derive for SDSS J0815+4304 
(which is a good Compton-thick candidate) is rather 
modest i n comparison with values > 1 seen in some 
ULIRGs ()Spoon et all l2002[ ) which are therefore likely 
to have even higher column densities. 

Water ice absorption and organic molecules (including 
the stretching modes of CC bonds in ring-like aromatic 
components) can contribute at wavelengths > 6/im. The 
precise shape of the opacity curve at 6 — 8/im is difhcult 
to determine because it is affected by the choice of the 
continuum and by assumptions about the strength of the 
7.7/im PAH feature. Assuming a power-law continuum 
and no PAH emission, we derive the opacity curve shown 
in Figure [5)3. For the organic component, a broad range 
of opacity curves h as been obtained in the laborato ry at 
these wavelengths (jPendleton fc Allamandolall2002[ ). but 
none appears to have the required shape and strength rel- 
ative to the narrow aliphatic features (a typical opacity 
curve of laboratory hydrocarbons is shown as the solid 
line in Figure [5)3) . A water ice feature'^ with a peak op- 
tical depth T[6/^m] = 0.30 can reproduce the absorption 
at A < 6/im, but fails at longer wavelengths (dotted line 
in Figure[5)3). More importantly, there is no evidence for 
the stronger 13/im water ice absorption feature in our 
spectrum (T[13/im]/T[6Mni] > 2 for the water ice opac- 
ity curve; iGerakines et al.lll995f ). which rules out water 
ice as the dominant source of opacity at 6/Ltm. The spe- 
cific carrier of the excess 6.0 — T .S/xm opacity (whi ch has 
also been observed in ULIRGs. ISpoon et al.l[200l 120041 : 
iDartois et aIll2007D remains unidentified. 



3.4. PAH emission features 

Many Galactic and extragalactic objects show strong 
broad emission features at 3-13/im which are tho ught to 
originate fr om small carbonaceous dust particles (|Drainel 
I2OO2I [2OO3I : see top panels of Figure[T]for the template of 
this emission). In contrast to the aliphatic (chain- like) 
molecules responsible for the absorption bands discussed 
in the previous section, these emission features are due 
to polycyclic aromatic hydrocarbon (PAH) molecules in 
which < 50 carbon atoms are arranged on a hexago- 
nal lattice. The different feat ures arise from different 
C-C and C-H bending modes (jAllamandola et all 119891 : 
iDraind [20031 ). and the precise spectral positions and ra- 
tios of the bands depend on which specific molecules are 
present, on their ionization state and on the spectral 
shape of the illuminating emission. 

PAH emission is powered by UV photons, and there- 
fore can in principle be powered either by star-forming 
regions or by the active nucleus. Several arguments fa- 
vor the hypothesis that PAH emission predominantly 
traces star formation rather than AGN activity. In 

^ Water ice opacity taken from the database of 
the Sackler Laboratory for Astrophysic s in Leiden, 
|http: / / www.strw.leidenuiiiv.nl/~lab/databases/ 1 



particular, 7.7/im PAH emission is spatially extended 
in nearby AGNs and is suppressed near the nucleus 
()Le Floc'h et al.|[200lD . presumably because of the de- 
struction of PAH molec ules by ext reme-UV and X-ray 
photons from the AGN (jVoitl 119921. Furthermore, the 
near-to-far-IR spectral energy distribution of AGNs with 
PAH emission is accurately represented by a sum of a 
non-P AH AGN templ ate and a star-burst galaxy tem- 
plate (jShi et all l20Cl7f l. Finally, for quasars with CO 
emission, the CO-PAH correl ation is in ag reement with 
that of star-forming galaxies (jShi et al.| [2007'). 

Since PAH[7.7/im] is blended with PAH[8.6/im], Si ab- 
sorption, the [NeVI]7.63/im emission line and potentially 
ice/hydrocarbon absorption as well, we use the strength 
of the 11.3/im PAH feature as a star formation indicator. 
Based on the fits described in Section 13.11 six objects 
were judged to have a detected 11.3/im PAH feature. In 
Table [U we list the average of the fluxes obtained by the 
two methods (global fitting with PAHFIT and local fit- 
ting of a continuum-f Drude profile model) and assume 
that they are accurate to no better than 30%. We also 
assume that the PAH emission is not affected by Si ab- 
sorption (i.e., that it originates from outside the region 
affected by absorption). If this is not the case, then the 
11.3/i m lumin osities listed in Table[T]are underestimated. 

Lut z et al] ([1998) proposed a classification scheme in 
which they used the ratio of peak (continuum sub- 
tracted) PAH[7.7/zm] flux density to 7.7/im continuum 
to distinguish between starburst- and AGN-powered 
sources, an d a similar i d ea is employed in the 'fork 
diagram' of ISpoon et"an (|2007D . A low PAH/mid-IR 
continuum ratio is characteristic of AGNs, because the 
spectral energy distributions of star-forming galaxies 
peak at much longer wavelengt hs than tho s e of A GNs 
(Figured]). By the criteria of iLutz et all ()1998f ). all 
our objects are AGN-dominated in the mid-IR. The 
PAH/continuum ratio can then serve as a crude estimate 
of the relative contribution of star formation in the host 
galaxy to the total IR energy budget. Figure [6] shows 
an anti-correlation between the [OIII]5007A/[OII]3727A 
line ratio and the PAH[11.3/im]/continuum[27.5/im] 
ratio, as well as between [OIII]5007A/H/3 and 
PAH[11.3/im]/continuum[27.5/im], which are detected 
with about 95% confidence. These anti-correlations sug- 
gest that for quasars [OIII]/[OII] and [0III]/H/3 de- 
crease as the star formation in the host increases, in 
agreement with class ical work on optical line diagnostics 
(|Baldwin et al.lll981h . 

We roughly estimate the total luminosity as- 
sociated with star formation using the rela- 
tion log(LsF[8 - lOOO/im]/lO^°L0) ^ 1 16 x 
log(L[PAH11.3Aim]/10i°LQ) + 2.38 (|Shi et all [200I . 
These values are listed in Table [1] In Figure [7l we show 
the comparison between the distribution of estimated 
star-formation luminosity in our objects (histograms) 
and the luminosity function of star formation in other 
AGNs and in field galaxies from a variety of samples 
as described in the figure caption. The median star- 
formation luminosity of our sample (5 x lO^^L©, which 
includes objects with upper limits on PAH[11.3//m]) is 
higher than that of field galaxies and of any other AGN 
sample. 

As high as the star formation luminosities are, in all 
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but one case they are less than hah of the total infrared 
luminosity obtained by directly integrating the observed 
IR spectral energy distribution. Thus most of the IR lu- 
minosity is due to the dust-reprocessed emission from the 
quasar (the one exceptio n is discussed below) . Models of 
AGN obscuration (e.g., iPier fc KroTi3ll992l ) show that 
dust-reprocessed emission is anisotropic even at mid-IR 
wavelengths, and the luminosities that we obtain in type 
2 quasars by integrating the IR emission, multiplying by 
47r and subtracting the star-formation contribution un- 
derestimate the quasar bolometric luminosity, by a fac- 
tor which depends on the viewing angle (see Figure 4 of 
Pier fc Krolik'l992D. When we average model fluxes from 
Pier fc Krolik (199^) over all angles such that the central 
source is not directly seen, we find that this factor ranges 
from 2.2 to 4.3. This implies that the true bolometric lu- 
minosities of the quasars in our sample are likely much 
larger than the star-formation luminosities of their hosts. 
Furthermore, in all cases the observed radio luminosity 
is strong er than w hat is expected from star formation 
alone (Helou et al.|[l985 ) , and most of the radio luminos- 
ity is due to the quasar, even though all the objects in 
our sample (except SDSS J0812-I-4018) are radio-quiet 
(based on their positi on in the Lradio — i[OIII] diagram; 
IZakamska et al.ll2004[ ). 

In the radio-loud type 2 quasar SDSS J0812+4018, the 
star-formation luminosity estimated from PAH[11.3/im] 
is similar to the total IR luminosity. In this object, the 
11.3/im feature appears anomalously strong compared to 
the 7.7^m one, and therefore the star-formation lumi- 
nosity is probably overestimated. Although it is usu- 
ally thought that high-lumi nosity AGNs destroy PAHs 
in their vicinity (|Voitl 11993 ). the processes that deter- 
mine PAH feature ratios are not well-understood, and 
it is possible that the emission from the quasar itself 
boosts the 11.3/im feat ure in t his ca se, as is seen in some 
low-luminosity AGNs (jSmith etTal] l2007h . An anoma- 
lously strong 11.3/im PAH feature is seen in some other 
extragalactic obje cts (e.g., in a buried AGN studied by 
iSpoon et al.l[2f)nl . 

3.5. Atomic emission lines 

We detect a number of ionized elements through their 
fine-structure emission lines (the lines that we have con- 
sidered are enumerated in the top panels of Figure [1]) . 
Just as certain optical line ratios are diagnostic of the 
physical conditions in the medium where they are pro- 
duced and of the intensity and spectrum of the ioniz- 
ing radiation field (jVeilleux and Osterbrocklll987f ). some 
mid-IR line ratios can also be used as such diagnostics 
(jSturm et al. 2002; Groves ct al. 2006). 

With ionization potentials of 126 eV and 97 eV, the 
NeVI and NeV species are produced exclusively by the 
active nucleus, whereas Nell (ionization potential 22 eV) 
can be excited both by the nucleus and by massive stars. 
Therefore, the [NeVI] / [Nell] and [NeV] /[Nell] fine ratios 
are diagnostic of the AGN contribution to the line emis- 
sion. In addition, for a fixed shape of the AGN spectrum, 
these line ratios increase with increasing ioiiizatio n pa- 
rameter (i.e., radiation densitv. [Groves et al.ll2006l ). We 
detect the [NeVI]7.63/im emission line in 11 out of our 12 
spectra. In objects with PAH emission, identification of 
[NeVI]7.63/im is somewhat complicated by the fact that 
it is right on top of the PAH[7.7/im] feature. The median 



[NeVI]7.63/im/[NeII]12.8/im ratio in our obj ects is 1.3, 
simila r to the sample of Seyfert galaxies bv lSturm et al.l 
(|2002D . where the median is 0.9. The [NeV]14.3/im emis- 
sion line is also detected in 6 out of the 7 objects in 
which the wavelength coverage of this line is available. 
The median [NeV]14.3Mm/[NeII]12. 8/im ratio is about 
1, similar to that of the sample of ISturm et al.l ()2002D 
and somewhat larger t han that of Seyfert galaxies from 
iTommasin et al.l ()2008[ ) who found a median ratio of 0.5. 

The simple 'linear mixing' mod el of an AGN-f starburst 
spectrum bv ISturm et all (|2002D shows that our median 
value of [NeV]14.3/im/[NeII]12.8/im= 1 roughly corre- 
sponds to a 100% AGN contribution, and our lowest 
value of this ratio (0.5 in SDSS J0050-0039) imphes an 
AGN contribution of 50%. Since [NeV] and [NeVI] are 
produced in close proximity to the nucleus, they are more 
affected by extinction in the host galaxy or in the cir- 
cumnuclear material than [Nell]. Thus, our measured 
[NeV] /[Nell] and [NeVI] / [Nell] ratios are in fact lower 
limits, and so is the AGN contribution calculated from 
the 'linear mixing' model. We do not find a strong corre- 
lation between the Ne line ratios and the PAH equivalent 
widths. One likely reason is that a range of ionization 
parameters is present within our sample, so the depen- 
dence of Ne line ratios on the star formation contribution 
is washed out. 

The one emission line that shows significantly differ- 
ent properties in type 2 quasars and Seyfert galaxies 
is [SIV]10.51/im (ionization potential 35 eV). This line 
is strongly detected in 11 of our 12 objects, and the 
median ratio is [SIV]10.51/fm/[NeII]12.8/fm~ 1.2 (un- 
corrected for Si absorption). This ratio is significantly 
higher than that of the Seyfert galaxies in the sample of 
iTommasin et al.l (120081) ( media n:^ 0.4) and higher still 
than that of lSturm et al.l ()2002f ) (median~ 0.1). Perhaps 
this difference is due to the larger ionization parameter, 
since the quasars in our sample are significantly more 
luminous (median L(IR)= 10^^'^erg/s) than the Seyfert 
galaxies in the com parison samples (e.g., lO^'* '^ erg/s in 
ISturm et al.l[200l . 

3.6. Molecular hydrogen emission lines 

Rotational transitions of molecular hydrogen (H2) al- 
low us to directly observe warm (> 200K) molecular gas 
and esti mate its excitation tem perature and mass (e.g., 
Rigo poulou et a l. '2002"; ' Armus et aLH2006l : iHigdon et all 
i2006t iRoussefet al.. ,2007f r The molecule is symmetric, 
so the transitions are quadrupolar and therefore weak, 
and their detectability is limited by the low signal-to- 
noise ratio of our spectra. Furthermore, some of the hy- 
drogen lines blend with other spectral features, which is 
especially problematic because of our poor spectral res- 
olution. 

Despite these difficulties, we detect several rotational 
H2 lines in two objects (Tableland FigurejHl). The confi- 
dence of the detections, uncertainties in fluxes and upper 
limits are all calculated using simulations of weak unre- 
solved emission lines in noisy spectra. For a line S(J), 
the transition is from the upper level with angular mo- 
mentum J-l-2 to the lower level with angular momentum 
J. Under the assumption that emission is purely sponta- 
neous and without correcting for extinction, the line fiux 
Fj can be directly related to the population of the J + 2 
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level: 
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where Dl is the luminosity dis tance to the sourc e, A 
are the Einstein coefhcients (e.g.. lTurner et allllQ??! ) and 
— Ej = hc/X is the energy of the transition. The 
rotational energy levels are given by 

Ej = 85.35K•A:iJJ(J+l)-0.068K•fcBJ2(J+l)^ (2) 

where gj is the degeneracy of the levels, 2J + 1 for 
even J and 3(2 J + 1) for odd J (assuming equilibrium 
ortho-to-pa ra ratio) and /cb is the Boltzmann constant 
(jHuber fc Hcrzberg 197l|). In this expression, the first 
term is the quantum-mechanical analog of the energy of 
a rotator with angular momentum J. As J increases, 
the energy of the ground vibrational level is changed by 
rotational stretching, which is taken into account by the 
second term (Landau & Lifshitz 1991). If different levels 
are populated in a Boltzmann equilibrium, the line ra- 
tios are indicative of the excitation temperature which is 
defined as 

cxexp -— — . (3) 

gj \ kb^cxcJ 

In the spectrum of SDSS J1641-H3858, we detect 
S(3)9.665Mni, S(4)8.026Atm and possibly S(5)6.909Aim, 
allowing us to calculate an excitation temperature for 
these transitions. The S(l)17.03^m and S(4) hues 
are detected in SDSS J0050-0039 with a lower confi- 
dence. In Figure [HI we show the excitation diagrams 
for both objects, i.e., we plot ln{Nj/gj) from equa- 
tion II]) as a function of the energy of the J-th level. 
In SDSS J1641+3858, we find an excitation tempera- 
ture of about 700 K, somewh at higher than the tempera- 
tures found by Hi gdon et al. (2006) in a large sample of 
ULIRGs (300-400K). This difference should not be inter- 
preted as a difference in physical conditions, but rather 
as a selection effect - specifically, our lack of longer wave- 
length coverage precludes us from making observations 
of lower J transitions which trace colder gas. Indeed, 
in SDSS J0050— 0039 the excitation temperature derived 
from S(l) and S(4) is lower (~420K). 

Using these excitation temperatures we can calculate 
the mass of the gas. Levels are populated according to 
equation ^ with the normalization given by equation 
([T|). Summing up molecules with all J yields a mass of 
3.5 X IO^Mq in SDSS J1641+3858 (emitting at temper- 
ature 700K) and of 1.8 x IO^Mq in SDSS J0050-0039 
(emitting at temperature 420K). For comparison, in sev- 
eral ULIRGs with five or more detected H2 rotational 
lines ()Armus et al.l 120061 : iHigdon et all 120061) . the gas 
is found at a range of excitation temperatures, with 
- IO^Mq of warm gas (T > lOOOK) and - lO^M© of 
cool gas (T < 300K). There could be large amounts of 
even colder gas which would be invisible in rotational 
transitions of H2. For example, if one of our objects con- 
tains lO^^M© of cold gas with T < 150K, the population 
of J 7^ levels is so low that even the lowest energy tran- 
sition S(0)28.22/^m (J = 2 ^ J = 0) is weh below the 
typical noise l evel o f our data, as well as of the data by 
IHigdon et all (f2006h . Such cold gas can be traced only 
through lower energy transitions of other molecules, such 
as CO(l-O). 



Although II2 rotational line emission is often associ- 
ated with star formation (a nd then it is strong ly cor- 
related with PAH emission, iRoussel et al] |2007|) . there 
is no evidence for PAH emission in SDSS J1641-I-3858. 
If H2 emission originates in the circumnuclear region 
(for example, d < lOOpc), then the column density of 
the molecular gas responsible for the observed emission 



is rather high, N[li2] > 5 x 10^^^^^) 



-2 cm-2 



and 



yet it maintains an excitation temperature of 700K. As 
in SDSS J1641-I-3858, H2 emission in lower-luminosity 
AGNs is often in excess of what is expected fr om star 
formation in the host galaxy (|Roussel et al.|[2007l ) . There 
is at the moment no conclusive model for H2 emission in 
AGNs. 

4. COMPARISON WITH IR SPECTRA OF OTHER TYPE 2 
QUASARS 

In this section we compare the IR properties of opti- 
cally selected type 2 quasars with those selected by other 
methods. We discuss IR-selected (Section 14. 1|) . radio- 
loud (Section [42|) and X-ray selected type 2 quasars (Sec- 
tion [O]). 

4.1. IR-selected type 2 quasars 

Low-redshift (z < 1) quasars show roughly power- 
law spectral energy distributions in the IRAC bands, 
unlike nearby star-forming galaxies whose strong PAH 
emission dominates the 8/zm band. T he consequ e nt dif - 
fere nce in the inid-IR colors allowed iLacv et al] ()2004f ) 
and IStern et al.l ()2005f ) to select tens of candidate type 
2 quasars , which were then confirmed by optical spec- 
troscopy. iLacv et aTl (|2007t ) present IRS data for some 
of the objects from that sample; their Si absorption and 
IR colors are shown in Figure [2] as cyan squares. They 
also calculate the total PAH luminosity for the objects 
in their sample. On the basis of the PAH template de- 
scribed in Section 13.11 we estimate that the luminosity 
of the 11.3^m feature is 0.12 times that of the total PAH 
luminosity, and therefore the median 11.3/Km luminosity 
of the IR-selected objects is - the same as for 

our sample. 

Combining IRS and HST imaging dat a for the IR- 
selected type 2 quasars, iLacv et all (|2007[ ) find that the 
depth of the Si feature is correlated with the disk orienta- 
tion, in the sense that objects that appear more edge-on 
have a significantly deeper feature. These authors con- 
clude that some of the obscuration occurs in the plane 
of the host galaxy, rather than in the circumnuclear ma- 
terial on much smaller spatial scales (unless the circum- 
nuclear material is aligned with the plane of the host 
galaxy). In Figure [9] we show both samples (Lacy's and 
ours) on the 5*9.7 vs. axis ratio diagram. Morphologies 
and axis rati os are available from HST images for six o f 
our objects ([Zakamska et al.ir2006l : lArmus et all 119991 ). 
The IR-selected objects (open squares) tend to show that 
the Si absorption is stronger for the more inclined galax- 
ies, but the trend disappears when our optically-selected 
type 2 quasars are included (solid circles). There is a 
significant difference between the host properties of the 
IR-selected type 2 quasars and the objects in our sample. 
The former lie in disk-dominated, disturbed galaxies with 
dust lanes and active star formation. On the contrary, 
all our objects are hosted by bulge-dominated galaxies 
(2/3 of the hosts are ellipticals). Since the IR-selected 
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sample is hosted by disk galaxies, the axis ratio is a good 
measure of inclination in the IR-selected sample, but not 
in ours, and the lack of correlation in Figure [5] is unsur- 
prising. Additionally, [Nell] 1 2 . 8/im/ [Oil] 3727A ratio is 
on av erage 0.8 dex higher in the sample of iLacv et al.l 
()2007t l than in our objects. This difference is suggestive 
of a higher optical extinction in the host galaxies of the 
IR-selected objects than in ours, as supported by obser- 
vations of their optical morphology. 

Why is there such a difference between the host prop- 
erties of the two samples? The redshift ranges of both 
samples are similar, and so are the values of Si absorp- 
tion, IR colors and star-formation luminosities. Opti- 
cal properties (such as line ratios) are also indistinguish- 
able. The only significant difference seems to be in the 
luminosities of the two samples. Indeed, the median 
[OIII] 5007A luminosity of the objects from iLacv et aTl 
(|2007t) is log(L[OIII]/LQ) = 8.75 (calculated from their 
pubhshed [OII]3727 A luminosities using t he [Oil] vs. 
[OIII] correlation by IZakamska et al.l l2003f ) , while it is 
log(L[OIII]/io) = 9.41 for the objects in our sample. 
Similarly, the median 27.5/im monochromatic IR lumi- 
nosity for Lacy's sample is lO''^'^ erg/s, while it is 10^^'^ 
erg/s for our sample. It appears that there is a qual- 
itative change in the host properties of type 2 quasars 
occurring at a luminosity bracketed by the two samples. 
As it is possible that some subtle selection effect leads to 
the observed difference in host properties, this luminosity 
dependence clearly warrants further study. 

At higher redshifts (z > 1), type 2 quasars may be 
selected via a combination of IR and radio properties. 
Strong radio emission ensures that the object is an ac- 
tive nucleus (in the 'radio intermediate' regime) rather 
than a star- forming galaxy, while a large mid-IR/near- 
IR ratio is used to pre-select obs cured ra ther than unob- 
scured objects fjMartfnez-Sansigre et al.l[2006,) . The IRS 
spectra of ty pe 2 quasars selected using this method are 
presented bv lMartmez-Sansigre et al] ()2008[ ): about half 
of these objects were confirmed by optical spectroscopy. 
None of the objects show the 10/im Si feature in emis- 
sion, and roughly a third appear featureless, similar to 
our sample. Because of the higher redshift, the IRS spec- 
tra of these sources do not cover wavelengths longward of 
the 10/zm Si absorption feature, so we do not show these 
objects on the Si strength - IR color diagram, but the 
range of Si absorption in these objects seems to be com- 
parable to that of our s ample. The 'high-PAH' s u bsam - 
ple of the objects from iMartfnez-Sansigre et al.l ()2008l ) 
has PAH[11.3/Ltm] luminosities of 10^°~^°'^Lo - similar 
to the highest star-formation luminosities in our sample 
(Table [ij . This subsample appears to be well-matched 
in IR hmrinosity to the objects in our sample, while the 
' no-PAH' subsample is 3-4 times more luminous. 

iPoUetta et al.l (|2008l ) present a Spitzer spectroscopic 
study of very luminous {iyLi,[6^m\ > lO"*^'^ erg/s) high- 
redshift (z > 2) obscured quasars. These objects were 
selected from large-area Spitzer surveys based on their 
rest-frame NIR colors and a large IR-to-optical ratio, and 
therefore optical spectroscopy is not available for most 
sources. Those of the sources that showed lO/xm Si ab- 
sorption made up their sample of 21 obscured quasars. 
On the basis of this sample, the authors find that ob- 
scured quasars with Si absorption constitute more than 



half of the quasar population even at the highest lu- 
minosity end. Featureless IR spectra are common in 
the optically-selected and IR-selected samples of type 2 
quasars (and in the following sections we show that they 
are common in radio-selected and X-ray selected sam- 
ples as well), and therefore it is possible that yet more 
obscured quasars remain to be discovered at these high 
luminosities and redshifts. 

4.2. Radio galaxies 

In Figure O we show Si strengths and IR colors 
of rad io-loud AGNs from the sample of iHaas et aTl 
()2005f l. with filled magenta triangles for broad- 
line radio galaxies and filled blue triangles for 
narrow-line radio galaxies. Si strengths are esti- 
mated from the published spectra by computing 
ln(/y[9.7^m, observed] //i/ [9. 7/Ltm, extrapolated continuum]). 
Again, although a quarter of these objects show spectra 
with no Si absorption or emission, in the remaining 
cases the broad-line objects show Si in emission and the 
narrow-line objects show Si in absorption. The median 
[OIII]5007A luminosity of the narrow-line radio galaxies 
in this sample is log L [OIII] /L© — 8.56, so these objects 
are rough ly matched in lu minosity to the IR-selected 
sample of iLacv et al.l (|2007f) and are less luminous than 
ours. 

A larger sample of ra dio-loud AGNs with IRS obser- 
vations was compiled bv lClearv et al.l ()2007f ). More than 
half of their sample of 28 sources have no detectable Si 
feature. In the remaining objects. Si appears in emis- 
sion in radio quasars and in weak absorption in radio 
galaxies (which in the optical either show narrow lines or 
no lines at all). There are two exceptions to this trend: 
3C343, a radio quasar, shows Si weakly in absorption, 
and 3C325, a radio galaxy, shows Si weakly in emission. 
iGrimes et al.l ([^005) find weak broad lines and a red con- 
tinuum in the optical spectrum of 3C325, and they there- 
fore classify this object as a reddened quasar rather than 
a radio galaxy. Since the rest-frame UV opacity of dust is 
much larger than its IR opacity, a small amount of red- 
dening could produce strong UV extinction, while the 
IR spectrum would be transmitted unaffected. We fur- 
ther comment on the correspondence between the optical 
classification and the appearance of the Si feature in the 
next section. 

The range of Si absorption appears to be limited in 
both samples of radio-loud objects, with more than 
half showing no feature at all and no objects having 
5*9.7 ^ ^1- The distributions of 5*9. 7 in radio -quiet type 
2 quasars (taking our sample and that of Lacv et al. 
|2007) and in radio galaxies in the sample of iHaas et al. 

are different at > 90% level (determined using 
the Kolmogorov-Smirnov test). One possibility is that 
the contribution of jet to the mid-IR emission dilutes 
the spectrum, making the Si feature weaker relative to 
the total continuum. A larger sample of radio galaxies 
with well-measured Si strength is needed to investigate 
this difference in more detail. 

F ew, if any, obj ects in the samples of lHaas et al.l (|2005l ) 
and lClearv et al.l (|2007D show PAH emission . In a large 
samp le of 3CR objects (incl u ding those fromlHaas et al.l 
[2OO5I and lClearv et al.ll200l . lShi et all ([2007h find a me- 
dian star-formation luminosity of 3 x 10^^ Lq among ob- 
jects with detectable PAHs, and 6 x IO^Lq when the 
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upper limits on PAHs are accounted for. This level of 
star formation is lower than in PG quasars or 2MASS 
AGNs, and more than an order of magnitude lower than 
that of optically and IR-selected type 2 quasars. 

4.3. X-ray selected type 2 quasars 

IR spectra of seven X-ray sele cted type 2 AGNs were 
presented bv lSturrn et al.l (|2006f ). A wide ranee of IR lu- 
minosities is represented, from vL^lTfim] = lO'*'^'^ erg/s 
to lO'''^ '* erg/s. Unlike our sample, the IR spectra of 
these objects are mostly featureless. PAH features are 
firmly detected in only one object, and in all sources the 
star-formation luminosity is constrained to be < a few 
times 10^"Lq - significantly smaller than that of the ob- 
jects in our sample. None shows the 10/zm Si feature in 
absorption (in one of the seven cases the spectrum does 
not cover the wavelength range affected by the Si fea- 
ture). One object (LH 28B) may have Si in emission. 
This source shows broad lines with a large Balmer de cre- 
ment in its optical spectrum (jLehmann et al.ll2000[ ). so 
it may be another case of a reddened AGN, similar to 
3C325 described above. 

One possible reason for the difference in the appear- 
ance of the Si feature is that the two samples probe 
different ranges of column densities. Of the seven ob- 
jects in the X-ray selected sample, only one (CDF-S 202) 
has Nh > 10^'' cm~^, and that is the object in which 
the wavelength coverage of the Si feature is not avail- 
able. The remaining objects have column densities below 
]^q23.5 qyr~'^. Although we do not have direct measure- 
ments of the column density in our sample, more than 
half of our objects (and in particular, two objects with 
the deepest Si absorption) are probably Compton-thick, 
as we discussed in Section IX2l and as w e show in Figure H 
(jPtak et al.l 120061 : IVignah et aDl2006l ). None of the four 
type 2 quasars in our sample with Nh ^ lO^'^'^ cm~^ 
sho ws strong Si ab sorption (top right corner of Figure 
11]). IShi et all ()2006f ) also find that deep Si absorption is 
more likely to occur in Compton-thick sources. There- 
fore, a bias against Compton-thick sources in X-ray se- 
lected samples may result in a bias against deep Si ab- 
sorption in their IR spectra. 

5. SI STRENGTH VS. IR COLOR DIAGRAM 

5.1. Quasar sequence 

In this section we discuss the position of various types 
of sources (AGNs and ULIRGs) on the Sg.r-IR color dia- 
gram (Figured]). We start with the type 1/type 2 quasar 
sequence. In our HST data, some type 2 quasars clearly 
show conically shaped scattering regions in blue contin- 
uum light, with opening angles (projected on the sky) 
ranging from 5° to 60° (jZakamska et al.ll2006[ ). If we 
could determine the inclination of the scattering cones 
relative to the line of sight, we would know our viewing 
angle relative to the o bscuring material. For example, 
IZakamska eFal] ()2006f ) argued that in SDSS J1039+6430 
the scattering region to the south of the nucleus is ori- 
ented toward the observer because the opposing northern 
scattering region is fainter. If so, we may be able to see 
down to hotter areas where dust is heated by the quasar, 
even if not deep enough to see the broad-line region it- 
self. This would be consistent with the absence of a Si 
absorption feature in this object (15*9.71 < 0.05). On the 



other hand, the two objects that have rather symmetric 
scattering regions in the HST images, SDSS J1323— 0159 
and SDSS J1106-I-0357, may be seen almost exactly 
edge-on. Both these objects show the Si feature in ab- 
sorption (5g.7 = —0.16 and —0.23, correspondingly). 
SDSS J0920-I-4531, the only object in our sample with a 
possible Si emission, is also the only object in which an 
extended scatteri ng region could not b e identified from 
the HST images (jZakamska et al.ll2006D . It is likely that 
the compact blue source in the center of this object is 
the one responsible for the observed polarization, and 
it is possible that the scattering cone appears compact 
because it is viewed close to its axis (almost face-on). 
Thus, we may be seeing a dependence of Si strength on 
the viewing angle within the type 2 quasar population. 
We also saw in the previous section that the appearance 
of the 10/im Si feature (in emission or in absorption) is 
well-correlated with the optical classification (type 1 or 
type 2). These findings imply that the type 1/type 2 
quasar sequence in Figure [2] is driven at least in part by 
the viewing angle. 

This strong relation between the appearance of the 
Si feature and the viewing angle is rather surprising. 
First, the broad-line region is much more compact 

0.05pc • ( Sms-" )"' from the virial argument 

for bulk motio ns; cf. reverberation mapping sizes by 
iPetersoiil 1199,^ than the hot dust which is on scales 
larger than the region of dust sublimation 0.4pc • 

( io->6c^gs-i )^^^ (l^W) fro™ heating vs. cooling 
energy balance for dust particles] . Therefore, there might 
be plenty of lines of sight for which the broad line region 
is obscured while the hot dust region is not, with the Si 
feature appearing in emission, as s een in some geometrie s 
in radiative transfer models by iPier fc KrolikI (|1992l ). 
However, very few, if any, of the optically-classified type 
2 quasars show the feature in emission. The second com- 
plication is that the appearance of the Si feature is very 
sensitive to the geometry of the torus, its total optical 
depth and the distribution of dust within it. For ex- 
ample, models of Co mpton-thick tori alw ays show the 
feature in absorption (jPier fc KrolikllT992t l. independent 
of orientation, because the torus is too optically thick to 
be heated through in the direction parallel to the axis of 
symmetry, and large temperature gradients are present 

along all lines of si ght. 

In clumpy tori (jNenkova et al.l |2002[ ) the Si appear- 
ance depends mostly on the distribution and sizes of 
clumps, rather than on the orientation of the observer* 
. When dumpiness is reduced and the volume-filling fac- 
tor of clumps is increased. Si shows a stronger depen- 
dence on the observer's orientation (jSchartmann et aTl 
l2008f l. Therefore, the observed relation between the op- 
tical classification and the appearance of the Si feature 
may prove useful in constraining clumpy models of AGN 
obscuration. Anothe r characteristic of clumpy models of 
iNenkova et al.l (|2002f ) is that the shape of the Si feature in 
the final spectrum is often significantly different from the 
shape of the feature in the opacity curve, due to a super- 
position of emission from regions with different temper- 
atures. This diversity is not seen in those of our spectra 

* The library of spectral energy distributions is provided by the 
CLUMPY group at https://newton.pa.uky.edu/~clumpyweb/. 
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that show Si in absorption, while the very low color tem- 
perature of the possible Si emission in SDSS J0920-I-4531 
is a challenge to all existing models. The high fraction 
of featureless spectra, the centroid and the shape of the 
IG/im feature, the relative strengths of the IG/im and the 
18/zm features (not available f or most our objects , but 
employed in a ULIRG study of iSirockv et"all[2008h and 
the relationship bet ween Si strength and the observed 
Nh (jShi et al.ll2006[ ) may prove to be useful diagnostics 
of AGN models. 

The observed mid-IR F^[lA.h]/ F^[27 .b] color of aU 
quasars in Figure [2] has a median of 0.5±C).2 (with 60% 
of objects falling within this range). This is not well 
reproduced by either smooth or clumpy models, which 
typically predict i^^ [14.5]/i^y[27.5] > 1. However, the 
models do predict that this color should be rather sim- 
ilar in type 1 and type 2 quasars, as observed (type 1 
quasars in Figure [5] have a ratio of O.G^q'j, while type 2 
quasars have a ratio of 0.4 ± 0.1). In order to bring the 
model color in agreement with observations, one could 
postulate a cold dust component in the host galaxy. As 
long as the temperature of this component is < 80K, its 
flux density at 14.5/im is more than an order of mag- 
nitude smaller than that at 27.5/^m (and it makes an 
even smaller contribution to the wavelengths of the Si 
feature). We can estimate the minimal mass of this com- 
ponent by assuming that it is an optically thin emitter 
[Fy(A) oc i?i,(T, A)t(A)] and by requiring that it boosts 
the luminosity at 27.5/im (taken to be 10"*^ ei'g/s) by 
a factor of two (see Figure [2]), which is sufficient to 
bring the [14.5] /f^ [27 5] from 1 to 0.5. The Milky 
Way opacity curve from iDraind (|2002f ) is used to con- 
vert between optical depth and hydrogen column den- 
sity {Nh ^ 1.2 X 102^r[27.5/^m] cm'^). The minimal 
hydrogen mass is 3 x 10 M© if the cold component has 
T = 80K, but the required mass grows rapidly as the 
assumed temperature is decreased (e.g., 1.1 x lO^^M© 
at T = 50K), since more mass is required to produce a 
given amount of flux at a colder temperature. The mass 
may be underestimated by another factor of 4 if the gas- 
to-dust ratio is simila r to the SMC value rather than to 
the Milky Way value ()Draindl200§ ). 

Because the radiative transfer models of circumnuclear 
obscuration in AGNs predict similar colors for type 1 and 
type 2 objects, it is impossible to distinguish between 
multi-temperature circumnuclear obscuration and inter- 
vening cold absorption on the basis of IR spectra alone. 
Indeed, the opacity curve of dust is rather wavelength- 
independent (outside the strong Si features at 10 and 
ISfim), so intervening cold absorption does not change 
the colors in a significant way (see the 'screen of cold 
dust' absorption vector in Figure [2^). In other words, 
the IR spectra of type 2 quasars look just like they would 
if there was a screen of cold dust in front of a type 
1 quasar. The correlations between S'9.7 and IR colors 
found in our sample of 12 objects shown in Figures [21d,c 
also seem consistent with the 'screen of cold dust' model. 
Our HST observations of scattering co nes strongly sug- 
gest obscuration close to the nucleus (jZakamska et al.l 
|2006[) . so perhaps the outer parts of the obscuring torus 
are acting as this screen, while contributing emission at 
long wavelengths as described in the previous paragraph. 
Obscuration by intervening cold dust in the host galaxy 



seems to be imp ortant in the less luminous sample of 
iLacv et aD (|2007( ). as we described above. 

In type 1 quasars, there is a well-defined maximum 
value of the strength of Si emission, around S'9.7 = 0-5. 
To understand the origin of this maximum, we use 
DUSTY'"^ to calculate emission of optically-thick slabs 
illuminated by the quasar spectrum on one side. All 
incident radiation is absorbed in a thin layer near the 
illuminating source. Some of it is transmitted through 
the slab, while some of it is re-emitted toward the source. 
The spectrum of this 'reflected' component is almost in- 
dependent of the optical depth, but depends somewhat 
on the assumed temperature at the illuminated surface, 
with Si strength S'9.7 = 0.9 - 1.6 x log(ri„/1200K). So 
if the temperature on the illuminated surface is close to 
the dust sublimation temperature (~ 1200K), the model 
reproduces the observed maximum Si strength. The ob- 
served emission from quasars includes emission from the 
much colder outer parts of the obscuring material and 
thus is colder than the predicted reflected emission from 
the slab {F^[U.5]/ F^[27.5] ~ 2 at Tin = 1200K). 

5.2. Seyfert galaxies and ULIRGs 

The Seyfert galaxies in the sample of lHao et al.l ()2007D 
show a range of Si absorption and emission similar to 
that of the higher-luminosity quasars, as well as similar 
IR colors (F^ [14.5] /F^ [27.5] = 0.44tg;^^). However, the 
optical classification of Seyferts does not correlate with 
the appearance of the Si feature nearly as well as it does 
for quasars. In particular, more than 50% of Seyfert 1 
galaxies show Si in absorption, while only 10% of type 
1 quasars do. One possible explanation for this differ- 
ence is that even a modest amount of star formation in a 
Seyfert galaxy can significantly modify its IR spectrum, 
since the nucleus is not all that luminous. This possibil- 
ity is illustrated with a dotted line in Figure [21 where we 
show a locus of linear combinations of quasar and ULIRG 
spectra to mimic an AGN in a strongly star-forming 
galaxy. We see that when a star formation contribution 
is increased, the IR colors rapidly become much colder, 
contrary to what is ob served. Indeed, Seyfert 2 galax- 
ies from the sample of iHao et al.l (|2007[ ) have a median 
color Fj.[14.5]/Fj.[27.5] = 0.38 and a median Si strength 
S9.7 = —0.41, so they are warmer than the composite 
quasar -t-ULIRG objects on the dotted curve. Another 
possible explanation for the difference between IR prop- 
erties of Seyfert galaxies and quasars is that there is some 
systematic geometric or physical difference between the 
two classes because of their different bolometric luminos- 
ity. For example, it is possible that the opening angle of 
obscuration or the degree of dumpiness depend on lumi- 
nosity, resulting in different relationships between the IR 
properties and the viewing angle in Seyferts and quasars. 
Further theoretical work may clarify this interesting dif- 
ference. 

Type 2 quasars are significantly warmer than are 
ULIRGs for a given strength of Si absorption. Radia- 
tive transfer models show that this difference is not due 
to the different shapes of the illuminating spectra, since 
all incident em ission is absorbed with in a thin layer close 
to the source (jLevenson et al.l 120071 ). The total optical 

^ The DUSTY code is available at 
|http://www. pa.uky.edu/~moshe/dusty/ 1 
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depths are also thought to be comparable in quasars and 
ULIRGs. Instead, the differing IR colors probably reflect 
different geometric distributions of dust in quasars and 
ULIRGs. If the energy sources in ULIRGs (starbursts 
or AGNs) are completely embedded (that is, covered by 
obscuring material along all lines of sight), we would not 
see any warm dust. Using DUSTY, we explore optically 
thick spherical shells illuminated on the inside for a range 
of shell geome tries, optical depths and mass distribu- 
tion, following iLevenson et al.l (|2007f ). We find that the 
strength of Si absorption seen in ULIRGs is reproduced 
by models with IR optical depth, T[9.7^m], of order a 
few. There is a trend of colder IR colors for less cen- 
trally concentrated shells, but all model colors (for rea- 
sonable Si strengths) are warm, with F^, [14.5] [27.5] 
between 0.1 and 1. It appears that a significant fraction 
of ULIRGs (those with colder colors) require a cold dust 
component just as AGNs do. 

6. CONCLUSIONS 

We have conducted Spitzer IRS spectroscopy of ten 
type 2 quasars and obtained spectra of two more objects 
from the Spitzer archives. The sample was selected from 
the SDSS spectroscopic database on the basis of optical 
emission line properties and was required to have a min- 
imal [OIII]5007A luminosity of lO^L© and a minimal IR 
flux to obtain good quality spectra, but was otherwise 
unbiased with regard to any multi-wavelength proper- 
ties. In particular, the sample appears to represent a 
wide range of column densities along the line of sight, 
and half of the objects are likely to be Compton-thick. 
Optically-selected type 2 quasars show a diversity of IR 
spectral properties. Most objects show high-ionization 
fine structure emission lines such as [NeVI]7.63/im, essen- 
tially ruling out star formation as the dominant energy 
source. The lOfim Si feature appears in absorption in 
eight of the objects with a large range of strengths, and 
in another three the feature is not detected. In the one 
remaining object (SDSS J0920-I-4531) the feature may 
have been detected in emission, peaking at 11/im, which 
corresponds to an unusually low color temperature of 
2± llOK. 

Half of the sample show PAH features. One type 2 
quasar (SDSS J0815-I-4304) shows absorption due to or- 
ganic carbonaceous grains and possibly water ice absorp- 
tion. We detect three rotational lines of molecular hydro- 
gen in one object (SDSS J164H-3858), with excitation 
temperature Tcxc = 700K and mass M = 3.5 x lO'^M©. 
Molecular hydrogen is also detected (with a lower confi- 
dence) in SDSS J0050-0039. 

The median IR star-formation luminosity of type 2 
quasar host galaxies inferred from PAH emission is 

5 X lO^^L©, but the bolometric output of our objects is 
dominated by the obscured quasars. These levels of star 
formation are similar to those found bv lLacv et al.l ()2007l ) 
and iMartfnez-Sansigre et al.l ()2008f ) in other samples of 
type 2 quasars. The high star-formation luminosities of 
type 2 quasars - much higher than in field galaxies 
suggest that quasar activity is associated with intensified 
star formation in its host galaxy, as sugg ested by merger - 
driven models of quasar evolution (Hopk ins et al.|[2006f ). 
Moreover, the median star-formation luminosity that we 
find is higher than that of type 1 quasar host galaxies, 

6 X 10^" — 1.4 X IO^^Lq, depending on the quasar luminos- 



ity (|Shi et al.ll2007[ ). This trend supports the suggestion 
that the obscured fraction evolves with time, in the sense 
that the probability to see the q uasar as obscur ed de- 
creases with time after the merger (jHopkins et al.| [2006). 
In this case type 2 quasars would appear on average at 
an earlier stage of the merger and therefore with a higher 
star-formation rate than type 1 quasars, as observed. 

Objects in our sample occupy luminous bulge- 
dominated or elliptical hosts, whereas IR-selec ted type 
2 quasars from the sample of 'L acv et al.l (|2007f ) inhabit 
dusty disk galaxies. The two samples have different lu- 
minosities, with the objects from our sample being about 
0.5 dex more luminous {i'L^[27^m] = ip^S-S erg /s i n our 
sample vs. lO'^ '^ -^ erg /s in the sample of lLacv et al.l2 007'). 
iHopkins et all (pOOS) argue that at z < 1 most AGN ac- 
tivity at low luminosities is due to gas instabilities in the 
host galaxy, while high luminosity quasars are fueled in 
major mergers which result in formation of bulges. One 
possible explanation for the observed difference between 
host properties is that the two samples bracket a transi- 
tion between different quasar fueling mechanisms. 

Both optically-classified type 2 and type 1 quasars of- 
ten show spectra with no sign of the lO/xm Si feature. 
The fraction of such objects depends on how the sam- 
ple was selected; overall, about a third of all type 1 and 
type 2 quasars discussed in this paper do not show the 
feature. Most quasars that show Si in absorption are op- 
tically classified as type 2, while those with Si emission 
are type 1. Some models of clumpy AGN tori show that 
the appearance of the feature is primarily determined by 
the geometrical distribution of the material and its to- 
tal optical depth, while the viewing angle acts as a sec- 
ondary parameter. Therefore, the observed correlation 
between the optical classification and the appearance of 
the IR spectrum may be a useful discriminator of models 
of AGN obscuration. We further find that the correla- 
tion is significantly weaker in Seyfert galaxies than it is 
in quasars - specifically, half of Seyfert 1 galaxies show 
the lO^m Si feature in absorption. This suggests that ei- 
ther the geometry or the physical properties of obscuring 
material vary with luminosity. 
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TABLE 1 
Data summary 
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2 log 


, -L[OIII] 

' ^0 


S9.7 
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-0.47 
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46.4 


SDSS 


J005621. 72+003235. 8 
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45.25 
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40.9 
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11.9 


45.6 
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A A Q /1 ^ 
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44. / 
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44.00 


Ar\ Q 
4u.y 






AK 

40. u 


SDSS 


J081253.09+401859.9 


0.551 


9.39 


-0.14 


44.2 (1) 


44.99 


44.95 


43.2* 


9.7 


12.0* 


45.3 


SDSS 


J081507.42+430427.2 
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9.44 


-2.21 


<42.5 (3) 


45.37 


45.31 


40.8 


<9.2 


<11.5 


45.8 
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J091345.49+405628.2 


0.441 


10.33 


-0.32 


43.9 (1) 


46.48 


46.48 


40.8 
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<11.6 


46.7 


SDSS 


J092014.11+453157.3 
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9.15 
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45.13 
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0.00 
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45.78 
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40.2 
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46.3 
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J132323.33-015941.9 
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-0.31 
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39.8 
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<11.3 
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J164131. 73+385840.9 
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10.04 


-0.45 


44.8 (2,3) 


45.43 


45.46 


40.8 


<9.6 
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45.8 



Note. — All luminosities are K-corrected to the rest-frame. Lx = 47rZ)|^i<jf [2-lOkeV] (uncorrected for absorption); L[A] = inD^uFi^lX]. In the last column, 
L(IR) = 47rD|^F(IR), where i^(IR) is the integral of the IR spectral energy distribution from Spitzer photometry (Zakamska et al., in prep.). PAH luminosities are 
given uncorrected for absorption; star-formation luminosities are derived from PAH luminosities as described in Section l3.4l Typical uncertainties are 0.1 dex for IR, 
0.05 dex for radio and [OIIIJ5007A luminosities. Radio data are from FIRST survey (Becker et al. 1995^ Rest-frame 2—10 keV X-ray fluxes are obtained by fitting 
an absorbed power-law spectrum to the data. When only a few counts are detected, a P = 1.8 power-law spectrum is assumed. For Lx ^ 10*'^, Lx is determined 
to no better than 0.3 dex because of the Poisson noise and uncertainti es in the sp ectral shape. Refere nces fo r X-ray luminosities: (1) - unpublished (Ptak et al., in 
preparation, based on our recent Chandra and XMM programs); (2) - IPtak et al.l (2006); (3) - Vignal i et al.l (120061 ). 

SDSS J0812+4018 is the only radio-loud object; the rest are radio-quiet, given their position in the L[OIII]/Li.adio diagram dZakamska et al.ll200l) . The star- 
formation luminosity in this object is likely over-estimated, see Section |6]** SDSS J0920+4531 possibly shows the Si feature in emission, see Figure[3] If this is a 
correct interpretation, the peak strength of this feature is Sii^m = 0.13; however, there is no excess fiux from this emission at 9.7/xm, and therefore we list 59.7 = 0. 
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TABLE 2 

Molecular hydrogen lines 





S(l) flux, 


S(2) flux, 


S(3) flux. 


S(4) flux. 


S(5) flux, 


Object 


xlO"-'^ erg/s 


xlO~^^ erg/s 


xlO"-"-^ erg/s 


xlO^-"-^ erg/s 


xlO~^^ erg/s 


name 


(17.04^im) 


(12.28/im) 


(9.67/im) 


(S.OSAtm) 


(6.9lAtm) 


SDSS J0050-0039 


5±2 


<5 


<3 


1.4±0.7 




SDSS J1641+3858 


n/a 


<2.1 


1.0±0.5 


0.9±0.2 


0.6±0.3 



Note. — AU upper limits are at 95% confidence level (in the sense that 95% of lines with that 
flux would have been detected at 95% confidence level given the local signal-to- noise ratio). Spectral 
coverage of S{1) is not available for SDSS J1641+3858. We do not give upper limits for S(5) for SDSS 
J0050-0039 because it is blended with [ArII]6.985Atm, and in SDSS J1641+3858 it is only tentatively 
detected. There is also a possible detection of S(3) in SDSS J0920-I-4531 (Figure |3)l; in this case, the 
line falls in the high noise wavelength regime right between the SLl and LL2 modules of the IRS, so 
we do not provide quantitative measurements. 
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Fig. 1. — IRS spectra of type 2 quasars, plotted as vFy/c in unit s of mjy/^tm vs. re st-frame wavelength. The top- most panels show 
a PAH template derived from a composite spectrum of galaxies from [Smith et aTl II2007I) . the Si opacity curve from [Kemper ct al. (200'3) 
(dotted curve), and the positions of all the fine-structure atomic emission lines and rotational lines of molecular hydrogen (labeled SI, S2, 
etc.) that we have considered. Each spectrum includes object ID, redshift and the apparent optical depth of the 9.7/im Si feature derived as 
described in the text. Objects with detected PAH features are in the left column and those with no PAH features are in the right column. 
Strength of Si absorption increases from top to bottom. 
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Fig. 2. — Empirical relationship between the strength of the 10/.tm Si feature (59.7) and rest-frame IR colors for different types of objects. 
Black solid circles with error bars are SDSS type 2 quasars from this paper. There are four objects for which Sg.y = (the upper limits 
are all < 0.1). The vertical error bars are from comparison of measures of Sg.7 by different methods and the horizontal error bars assume a 
20% uncertainty in rest-frame colors (which includes the uncertainty in the photometric measurements and the uncertainty associated with 
the rest-frame correction). These error bars are representative of the measurements taken from other samples. Panels (b) and (c) show the 
IR colors which correlate most strongly with Sg.y for our sample (both correlations detected at > 98% significance). Additional notation 
pertains to Section [5] In panels (a)-(c) we show a 'screen of cold dust' absorption vector based on the extinction curve of .Kemper et aTl 
J2OO4). In panel (a), we show the cold dust (T < 80K) emission vector which corresponds to doubling the 27.5/.tm luminosity. The dotted 
line in panel (a) shows a locus of linear combinations of spectra cor r espon ding to the end points of each curve. For this particular realization, 
we take the c omposite spectrum of type 1 quasars from lHao et al.l l|2007l ) as the top right point and the spectrum of IRAS 01166-0844 from 
llmanishi et al.i t2007i^ as the bottom left point. 
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Fig. 3. — The spectrum of SDSS J0920+4531 showing possible Si emissi on, with PAH featu res (grey histogram) and without PAH features 
(black histogram; obtained by subtracting a scaled PAH template from [Smith et ani2007l) . The black point with error bars is based on 
our MIPS- 24 photometry of this source. A linear continuum (anchored at 8 /jm and 13;jm) has b een subtracted to show deviations more 
clearly. A similarly normalized spectrum of optically thin dust emission with lKemper et al.l I I2004I) opacity curve and temperature llOK is 
shown for comparison. Models with a lower temperature have a narrower Si feature and over-predict long-wavelength conti nuu m; models 
with a higher temperature ha ve a Si feature shifted to shorter wavelengths. Atomic emission features are indicated (Section 13.51 1: the 5(3) 
rotational line of H2 (Section 13.61 1 falls right between SLl and LL2 modules of IRS and is therefore only a tentative detection. 
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Fig. 4. — Strength of Si absorption plotted vs. the ratio of X-ray flux to that of [OIII]5007A. The 2-lOkeV X-ray fluxes have been 
corrected to t he rest-frame. The soh d line shows the mean ratio seen in Seyfert 1 galaxies and the dashed lines show the la range of this 
measurement l IHeckman et al. ]|200^). We then take a uFy =const X-ray spectrum and apply photoelectric absorption (S.Davis, private 
comm.) and Thompson scattering to calculate the X-ray/[OIII]5007A ratio in presence of line-of-sight gas with hydrogen column density 
Nji and solar metallicity (grey dotted lines); [OIIIJ5007A is assumed to be unaffected by absorption. 
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rest-frame A, fxm rest-frame A, /im 



Fig. 5. — Optical depth of the absorption complex in SDSS J0815+4304 (black histogram; grey histogram shows the subtracted 
[NeVI]7.63^m and [ArII]6.99^m emission lines). Left: a close-up view of the aliphatic features. Right: the opacity of the whole complex 
when the continuum is assumed to be a powe r-law betw een 5.5^m and 7.8/im. The thin solid line shows a typical laborato ry hy drogenated 
amorphous carbon opacity curve (copied from lDartois et al.l2007!) . The dotted line shows opacity of water ice at T < lOOK l lGerakines et ahl 
[19951) . ^ 
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Fig. 6. — Optical emission Unc ratios (open squares: [OIII]5007A/H/3; filled circles: [OIII]5007A/[OII]3727A) vs. the relative contribution 
of PAH[11.3/jm] to the IR emission. The latter is quantified by the ratio of PAH[11.3/.tm] luminosity to the monochromatic luminosity 
i^L^[27.5/im]. Both optical line ratios tend to decline as the relative PAH contribution increases - that is, as the relative contribution of 
star formation increases (anti-correlation is detected at the 95% confidence level using the rank correlation coefficient; correlation detected 
both when upper limits are included and when they are excluded). No definitive trends are seen for any other optical or mid-IR line ratios. 



20 




Fig. 7. — Distribution of star-formation luminosities in type 2 quasars (right vertical axis; solid line histogram — all objects; grey shaded 
histogram - upper limits) plotted with the same horizontal scale as the luminosity functions (left vertical axis) of s tar formation in local 
field galaxies (dotted line, ILe Floc'h et al. | 2005, arbitrarily normalized), in local Seyfert 1 galaxies (dashed line, IMaiolino et al.l 119951 . 
arbitrarily normalized), in Mg < —21.0 mag PG quasars (thin solid line. IShi et al.ll2007l : median star-formation luminosity 6 X IO^^Lq) 
and in Mg < —23.5 mag PG quasars (thick solid line, IShi et all 120071 ; median star-formation luminosity 1.4 X IO^^Lq). The median 
star-forma tion luminosity i n type 2 quasars (5 X lO^^L©) is higher than the median star-formation luminosity in any other group of AGNs 
studied by IShi et aEI | |2007|) . 
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Fig. 8. — Top row: rotational lines of molecular hydrogen in SDSS J1641-I-3858. (a-c): Solid lines show our best continuum -|- one 
Gaussian fits, or two Gaussians in the case of the S(5)-f [Aril] blend in panel (a), (d): Excitation diagram with the vertical axis reflecting 
the actual number of emitting molecules. The excitation temperature is by definition the inverse of the slope of the excitation diagram. 
Bottom row: tentative H2 line detections and excitation diagram for SDSS J0050— 0039. 
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Fig. 9. — Si strength vs. axis ratio of tlie liost galaxy. Solid circles are type 2 quasars from our sample. Axis ratios are taken from 
HST i mages bv lZakamska et al.l II2006I ) for five objects; for SDSS J0913 +4056 the ratio is estimated to be 0.5 from images by A rmus et al.l 
lfl999h . Open squares are IR-selected type 2 quasars from the sample bv lLacv et al.l I I2007I) who see a trend of increasing Si absorption with 
increasing e lliptic ity in their sample. When the objects from our sample are added, the trend disappears. The objects in the sample of 
ILacv et al.l l|2007l ) reside in dusty disks, whereas the objects in our sample are about 0.5 dex more luminous and reside in bulge-dominated 
galaxies. 



